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The cascade mass reconstruction approacii was applied to simulated production of the lightest 
stop quark at the LHC in the cascade decay p^tit— >X2tt^£_R^tt— )-Xi^^tt with top quarks 
decaying into hadrons. The stop quark mass was reconstructed assuming that the masses of gluino, 
slepton and of the two lightest neutralinos were reconstructed in advance. 

A data sample set for the SU3 model point containing 400k SUSY events was generated which cor- 
responded to an integrated luminosity of about 20 fb~^ at 14 TeV. These events were passed through 
the AcerDET detector simulator, which parametrized the response of a generic LHC detector. The 
mass of the ti was reconstructed with a precision of about 10%. 



I. INTRODUCTION 

If supersymmetry exists at an energy scale of 1 TeV, 
the study of third generation sleptons and squarks at the 
LHC is of a special interest. Their masses can be very 
different from that of sparticles of the first and second 
generation, because of the effects of large Yukawa and 
soft couplings as can be seen from the renormalization 
group equations. Furthermore they can show large mix- 
ing in pairs (fi, ii?), and (fL,f/j). A detailed 
discussion of possible SUSY effects at the LHC is given 
inQ. 

In this paper we consider the mass reconstruction of 
the lightest stop quark {ti) in the cascade decay 

g ^ tit ^ x"" ^ ^ii^tt ^ X?^^tt (1) 

with the top quarks decaying into hadrons. The gluino 
decay chain ([T]) is represented in Fig. [U in which all final 
state particles are explicitly shown. Here, the considered 
leptons are electrons and muons {i = e, The lightest 
neutralino x\ is invisible to the particle detector, whereas 
b quarks and light quarks (labeled as 'q' in Fig. [T]) are 
observed as jets. 
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FIG. 1. A gluino cascade decay chain with a stop quark pro- 
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Approaches to stop quark mass reconstruction in dif- 
ferent decay chains and for different points in the MSSM 
parameter space are discussed in literature (for example 
in 043)- Recent limits from searches for stop quarks 
were published in Q. 

The reconstruction of a SUSY event is complicated be- 
cause of the escaping neutralinos and of the many com- 
plex and competing decay modes. At present, there are 
two different approaches to SUSY mass reconstruction. 
The endpoint method, which has been widely studied 0- 
0, [9l-[l7j , looks for kinematic endpoints of invariant mass 
distributions. The second method is the mass relation 
approach [l8l - l2]| , based on the "mass relation equation" 
which relates the masses of the SUSY particles and the 
measured momenta of the detected particles. It was 
shown in that the mass relation approach can be 
successfully used for integrated luminosities as low as a 
few fb-i. 

In this work, the mass relation approach of [l^l is 
used for measuring the mass of the lightest stop quark 
tl, assuming about 20 fb~^ of integrated luminosity in 
LHC proton-proton collisions at ^/s = 14 TeV, under the 
assumption that the masses of the gluino, slepton and 
of the two lightest neutralinos have been reconstructed 
in advance with 10-20% uncertainty. At such low inte- 
grated luminosity, the stop mass reconstruction is quite 
challenging because of a high level of SUSY background 
(sparticles created in decay chains different from that of 
Fig.[T|) and Standard Model tt background. 

In this paper, a particular example is chosen (the SU3 
model; section |ll| to illustrate the method. First, an 
"event filter" is applied to suppress the background (sec- 
tion |lll|), making use of a likelihood function built with 
the known uncertainties on jet and lepton measurements, 
constrained by the mass relation equation. Next, events 
which do not satisfy the kinematic limits derived from 
the chain ([T]) are discarded (section ITV)) . Finally, a com- 
binatorial mass reconstruction method (section |V| is ap- 
plied to find the best estimate of the stop mass from the 
maximization of a combined likelihood function, which 
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TABLE I. mSUGRA parameters for tlie SU3 point [J. 



depends on all five sparticle masses (gluino, stop, slepton 
and the two lightest neutralinos) and is constructed for 
each possible permutation of the final state particles. 



II. SIMULATION 

The particular mSUGRA model chosen for this work is 
the bulk point SU3, one of the official benchmark points 
of the ATLAS collaboration 0, which is compatible with 
the recent precision WMAP data [22|- The values of 
the relevant mSUGRA parameters are given in Table ID 
For the SUSY particles in the cascade process ((J), the 
theoretical masses and the production cross section have 
been found by ISAJET 7.74 ^ as reported in Table HH 

The branching ratio for the gluino decay chain ([T]) at 
the SU3 point is 



- 25.2% 2 11.5% -0 11-4% 7 100% 

.9 — ^ ti — > X2 — > ^B. — > Xl 



0.33%. 



Monte Carlo simulations of SUSY production for the SU3 
model point were performed with the HERWIG 6.510 
event generator 



24 



Later, the produced events were 
passed through the AcerDET detector simulation [25| . 
which parametrized the response of a generic LHC detec- 
tor (ATLAS and CMS detector descriptions can be found 
in [2^ and (l^j). The efficiency for b-jet reconstruction 
and labeling was set to 80%, whereas the calorimeter re- 
sponse to electrons and jets was 

cr 0.12 



E 



VWgcv 
0.5 



E 



0.005 



0.03 



(2) 



(3) 



VWgcv 

For muons, the same response function as for electrons 
has been used as first approximation. 

A sample of 400k SUSY events (including signal and 
background processes) was generated. This approxi- 
mately corresponds to 20 fb""'^ of integrated luminosity 
for the SUSY SU3 point production cross section of 19 
pb at 14 TeV. The masses of g, X2; ^fli Xi listed in Ta- 
ble |IT] were given as input and it was assumed that these 
masses had been already measured with about 10-20% 
uncertainty with the method explained in (20j . 

In order to isolate the chain ([T]) and to suppress the 
background, the following selection cuts were applied to 
the reconstructed quantities: 

• exactly two isolated opposite-sign same-flavor 
(OSSF) Icptons (cither electrons or muons) with 
transverse momentum px > 20, 10 GeV; 



TABLE II. Theoretical masses and total production cross sec- 
tion a of SUSY particles at the SU3 point 



• two b-tagged jets with px > 50 GcV; 

• at least three jets with px larger than 150, 100, and 
50 GeV 

• at least nine jets with px > 10 GeV (including b- 
tagged jets); 

• no r-tagged jets; 

• Mcfi > 600 GeV and E'^^'^"" > 0.2Afoff, where M^e 
is the scalar sum of the missing transverse energy 
and the transverse momenta of the four hardest jets 
and E^^^^ is the missing transverse energy; 

• lepton invariant mass 50 GeV < Mgfi < 105 GeV. 

A total of 24 signal and 191 background events are left 
after applying these cuts: the SUSY background to the 
signal process ([T]) is thus significant (the classification of 
events as signal and background is based on the knowl- 
edge of the simulated information) . 

As shown in 0, the SM processes are suppressed sig- 
nificantly by the above requirements. The SM domi- 
nant background surviving the hard cuts is tt produc- 
tion, where both W bosons decay leptonically producing 
a hhll state. Since the 1 1 production cross section is 
about 833 pb, a 17M tt sample, corresponding to 20 
fb^"'^ of integrated luminosity, was generated with the 
HERWIG event generator. After applying the above se- 
lection cuts, only 21 1 1 background events survive. 

Note that the requirement of high hadronic activity is 
important for the 1 1 background suppression. If the cut 
on the total number of jets TVjet > 9 is loosened to 7 
jets, the number of 1 1 events surviving the selection cuts 
increases to 115. 

For every event, light jets (i.e. not tagged as b-jets) 
were combined in pairs whose invariant mass has been 
reconstructed. Only the independent pairs whose Minv 
is in the range 60-100 GeV (W boson region) have been 
retained and events with less than two jet pairs have been 
rejected. The invariant mass distribution obtained with 
all combinations is shown in Fig. [5J 

Next, all combinations of a jj pair and a b-jct have 
been considered, retaining only those with Minv =145- 
205 GeV, representing the acceptable W-b pairs in the 
top quark region (Fig. [J). After this step, from the initial 
215 SUSY events (signal -I- background), a total of 834 
b j j combinations were identified as candidates for a top 
quark decay t — >■ bW — )■ bjj, coming from 23 signal 
and 70 background events. Thus at this stage only one 
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FIG. 2. Light jet pair invariant mass distribution. Only jj 
pairs with mass in the range 60-100 GeV have been retained 
in the analysis. 

signal event was lost. On the other hand, applying the 
same procedure to the 21 tt background events left only 
3 events with 21 bjj combinations as candidates for a 
top quark decay. 

Naively, one would think that this selection should let 
all tt background events survive, but it is clearly not the 
case. If a W-boson decays leptonically an invariant mass 
of light jets is not related with the W-boson mass. Also 
a parton showering leads to redistribution of energy and 
this means that the momenta of the b-quark and the W- 
boson are not sufficient to give the top invariant mass. 
The final result is that most 1 1 events are dropped at this 
stage. 



III. BACKGROUND SUPPRESSION 

To suppress the background before the last (combina- 
torial) step, an event filter is used which assumes that the 
masses of ^Ri Xi s-re known (Table |ll]). The event fil- 
ter maximizes, for each event, a likelihood function con- 
strained by the mass relation equation, or equivalently 
minimizes the function: 

^ /^cvcnt ^mcas\2 

^^(m~,,mi)=Y,— f^^ + Xf{m,p) (4) 

where the index i runs over the two leptons and the two 
t-quarks, p™^'^^ and ai are the reconstructed momentum 
and its uncertainty, and is the true momentum, 

which is varied to find the minimum (only uncertainties 
in jet and lepton energy measurements are taken into 
account). The parameter A is a Lagrange multiplier and 
f{rh,p) = is the mass relation. 

The constraint f{m,p) = is the key of the "mass rela- 
tion approach" : it relates the masses of the SUSY parti- 
cles to the measured momenta of the final state particles 
in the chain ([T]) jT8l - l20j . The mass relation constraint 
is obtained as a solution of a system of four-momentum 
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FIG. 3. W-b invariant mass distribution. W-b pairs are se- 
lected in the range from 145 GeV to 205 GeV. 



constraints for each vertex containing SUSY particles in 
the decay chain (IT|). For example, for the gluino decay 
vertex one has to? = (p^o +ki-^ +ki2 +Pti +Pt2)'^ where the 
right hand part is in terms of the four-momenta of the 
lightest supersymmetric particle (LSP), the detectable 
leptons, and the reconstructed top quarks. Similar re- 
lations can be obtained for each SUSY vertex in the 
process ([1]), which contains four vertices, hence one gets 
four kinematic equations that can be solved to find the 
four components of p-^o in terms of the SUSY masses and 
the momenta of the detectable and reconstructed parti- 
cles. By substituting these components into the on-shell 
mass condition for the LSP, TO?n = pin , the mass rela- 
tion constraint includes all SUSY masses (to) and the 
momenta (p) of the detectable leptons and of the recon- 
structed t-quarks in the process ([T]). The explicit form of 
f{irL,p) = is given by [l8j-l20|. 

Note that in the decay chain ([Ij the locations of each 
of the two t-quarks and of each of the two leptons are un- 
known. Here, it is assumed that the t-quark with higher 
energy originates from the gluino decay. In addition, we 
assumed that the leptons with higher px originate from 
the X2 decay. The momentum resolution for a t-quark 
is computed as cr = crb-jot © fjot i © crjct 2, where all 
jets arc assumed to have the resolution mentioned above 
(equation ([3]) in section lll| . 

In the numerical minimization procedure by MC sam- 
pling, the gluino mass is left free to vary within ±20% 
from the value reported in Table |lT] (a Gaussian sam- 
pling), and the stop mass is left free to vary uniformly 
in the range 480±120 GeV because the stop mass is un- 
known at this step. The latter range has lower limit 
approximately related with the kinematic condition for 
a stop decay TO{^ > m^^o + nit and is limited from above 
by the sbottom quark mass, which is about 600 GeV for 
the SU3 point. 

For signal events, the event likelihood distribution has 
a maximum in the region of the {g, ti) mass plane cor- 
related with the true masses of g and ti. Hence, signal 
events should give a peak in the region of the true masses. 



4 




number of trials per event 

FIG. 4. The number of events versus the number of accepted 
trials per event with ^ < 10. The total number of trials is 
10^ 



On the other hand, for background events there is no 
strong correlation betvifcen the likelihood maximum and 
the true masses of g and ti. Therefore, if we arbitrar- 
ily chose a point in the (g, ti) mass plane in the range 
close to true masses, there is a very high probability that 
the value of ^ found from equation (jlj is smaller for a 
signal event than for a background one. For each event, 
10^ points in the mass plane were generated in the range 
720 ± 144 GeV for the gluino and 480 ± 120 GeV for 
the stop quark, and the was calculated. Results of 
event filter are presented in Fig. 01 This figure shows the 
number of events versus the number of trials per event 
in which ^ < 10. Two groups of events can be seen: a 
group with the number of accepted trials close to zero 
and a group with more than a few hundred of accepted 
trials. The first group presents background events for 
which unlike for signal events a minimum of need not 
be in the considered mass range. Thus if < 10 in at 
least 300 trials per event, this event was considered as 
a signal candidate and was retained for the subsequent 
analysis. 

Before the application of the event filter, we had 23 sig- 
nal events, 70 SUSY background events and 3 Standard 
Model 1 1 background events. After the application of the 
event filter, the SUSY background events were reduced 
approximately by a factor of 2 while a single signal event 
was lost: 22 signal events, 37 SUSY background events 
and 3 1 1 background events survived. 



type 

bB ^ bcC. (5) 

where particles A, B, C are invisible but particles b and 
c are considered as visible (they can be either directly 
detected or indirectly reconstructed from the properties 
of final state particles). 

In the literature, kinematic limits on the invariant mass 
distribution of 6 c pairs in decay ([5]) over a variable = 
{Pb + Pc)'^ often are given for the case in which at least 
one of the visible particles is massless. However, in some 
cases both particles b and c can have a non-negligible 
mass. For example, this is the case when a gluino decays 
into a stop quark and top quark. These kinematic limits 
for the case of all massive particles in process (O are 
given [H-d^ by 

q = \j{-R±^ -AQS)/2Q (6) 

with 

Q = Ml (7) 

R = (ml- Ml - Ml) {ml + m^) + 

+ {ml-M\ + Ml){Ml-ml^Ml) ^ ' 

S = Mliml^mlf + 

+ {Ml - Ml,) [ml {Ml ml - M^) + (9) 
+ ml {Ml + ml -Ml)] 

where the upper edge corresponds to the case when b 
and c particles are moving in opposite directions in the 
rest frame of particle A and the lower edge corresponds 
to the case when b and c particles are moving in the 
same direction. A nonzero lower limit is a consequence 
of nonzero masses of the particles. 

The kinematic limits from equation ([6]) for the process 
((H) are gmin = 375.1 GeV and (/max — 496.8 GeV for par- 
ticles created on-shell. Because SUSY particles can also 
be created off mass shell, we have set a wider range. By 
choosing a cut at 525 GeV in 1 1 invariant mass, the num- 
ber of events remaining after the event filter selection is 
reduced to 21 signal events, 34 SUSY background events 
and 3 1 1 background events. 

V. STOP QUARK MASS RECONSTRUCTION 



IV. KINEMATIC LIMITS 

It was shown 0-11, [i-[l3 that the endpoint method 
could be very useful in SUSY particle mass reconstruc- 
tion for finding relations between the masses of the SUSY 
particles involved in a decay chain and to determine their 
masses. Such a method allows mass reconstruction with- 
out relying on a specific SUSY model. In particular, the 
endpoint method can be applied to a decay chain of the 



If all sparticle masses but the stop mass in the decay 
chain ([T]) were known, the mass relation equation would 
allow finding the stop mass directly. However, we assume 
that these masses are known at integrated luminosity of 
20 fb~^ with an uncertainty of 10-20%. In order to take 
this into account, at the final step of the stop mass re- 
construction we allow for all sparticle masses to vary in 
ranges defined by their uncertainties. Because in this case 
for each event there are five unknown masses, at least five 
events are required to reconstruct sparticle masses. 
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FIG. 5. Reconstructed stop mass distribution including SUSY 
background and SM 1 1 background with integrated luminos- 
ity of 20 fb~\ The curve is the result of a Gaussian fit. 



The combinatorial procedure of [20| is used for the fi- 
nal stop mass reconstruction, applied only to the events 
that pass the event filter. It is important to make the 
most effective cuts in advance, because such procedure 
is computationally very intensive due to combinatorics: 
all possible groupings of five events in the sample are 
considered. 

For each set of five events, a combined likelihood func- 
tion is built and maximized. For sparticle masses the 
combined likelihood function for the combination is de- 
fined as the product of the maximum likelihood functions 
for individual events. Finding a maximum of the com- 
bined likelihood function for the combination is the same 
as searching for a minimum of the function 



i=l 



(10) 



In Eq. (Uni), min(x^ 



is a result of searching for 



a minimum of the X^ovcnt function for an individual 
event. For each of the five events in the combination, 
the min(x^ovont) fitted with 9 parameters (four particle 
momenta and five SUSY masses). 

The function for an individual event is defined by 



Xevent 



(P| 



p. 



mcas \ 2 



+ 



(11) 



mass widths are taken to be 15 GeV for the gluino, 2 
GcV for the stop and 1 GeV for all light masses. The 
first two numbers are comparable with the theoretical 
widths for the heavy SUSY particles. The last number 
takes into account the fact that light SUSY particles are 
either quite narrow or stable. We note that the results of 
the mass reconstruction are not strongly sensitive to the 
actual values of sparticle widths. In equation PT|) the 
mass relation and the il edge, which relates three light 
sparticle masses for chain ((T)), are included by means of 
the Lagrange multipliers Ai and A2. The tl edge (103.1 
GeV) can be obtained, for example, from Eq. ^ by plug- 
ging in the correspondent light sparticle masses and zero 
lepton masses. 

The minimization is done numerically by means of MC 
samplings of the parameter space. The sampling is uni- 
form for the stop mass, in the range: 480 ± 120 GeV. For 
the masses of g, Xj, Xi ^ Gaussian sampling is done, 
with mean values as given in the Table |IT] and standard 
deviations of 72, 20, 20, 20 GeV, which approximately 
corresponds to uncertainties of 10%, 10%, 15%, 20% in 
sparticle masses found in The MINUIT code [Hi 



is used to search for the minimum of the Xcomb func- 
tion (Uni). 

The reconstructed SUSY particle mass distribution, 
together with a Gaussian fit, is shown in Fig. [5] As can 
be seen in this figure, the reconstructed stop mass dis- 
tribution is described approximately by a Gaussian with 
the mean value of 483 GeV and standard deviation 36 
GcV. 

The asymmetry of the stop mass distribution toward 
higher masses can be explained by the nearness of the 
kinematic edge for the stop decay ti — >• X2 + t to the 
stop mass so that high masses are generated more often 
than lower values. This is not accounted for in the Monte 
Carlo sampling, for which we have uniformly generated 
the values of the stop mass. 

In order to understand the role of background and 
of the simulated detector effects, the stop mass has 
been reconstructed using generator-level momenta with- 
out SUSY background and Standard Model tt back- 
ground, and the result is shown in Fig. [51 In this case, 
as expected, the reconstructed stop mass is very close to 
the theoretical mass and the width of the distribution 
is smaller because of the absence of detector effects. It 
follows from comparison of Fig. [5] and Fig. [H] that the 
presence of background tends to shift the stop mass peak 
position to the region of higher masses. 



where the variables with superscript "event" are those 
with respect to which one has to minimize. 

The first term in equation ()lip takes into account the 
deviations of the reconstructed momenta of the t-quarks, 
or the measured momenta of the leptons, from the true 
ones. The second term takes into account the intrin- 
sic width of the SUSY particle masses, in the Gaus- 
sian approximation (instead of a Breit-Wigncr distribu- 
tion). The standard deviations corresponding to the 



VI. CONCLUSION 

We applied the cascade mass reconstruction approach 
developed in [2OI for reconstructing the mass of the 
generation supersymmetric quark, assuming 14 TeV 
proton-proton collisions at the LHC with integrated lumi- 
nosity of about 20 fb^^. At such relatively low integrated 
luminosity, the stop mass reconstruction is complicated 
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decay chain involving a stop quark. Our approach to the 
stop mass reconstruction is based on the consecutive use 
of an event filter and of a combinatorial mass reconstruc- 
tion method. 

In this work, wc considered the stop mass reconstruc- 
tion at the SU3 mSUGRA point and we obtained an esti- 
mate of the stop mass with a precision of about 10%. We 
expect that our approach should work for different MSSM 
parameters as well, provided that a decay chain contain- 
ing at least four successive two-body decays and involving 
five SUSY particles has a sufficiently large branching ra- 
tio to be identified in a heavy background environment. 



FIG. 6. Reconstructed stop mass distribution for true events 
without background with integrated luminosity of 20 fb^^. 
The curve is the result of a Gaussian fit. 



because of a high level of SUSY and Standard Model 1 1 
backgrounds and of the low branching ratio for the gluino 
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